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Abstract 

We propose a theory of two-magnon resonant Raman scattering from anti- 
ferromagnetic insulators, which contains information both on the magnetism 
and the carrier properties in the lighly doped phases. We argue that the con- 
ventional theory does not work in the resonant regime, in which the energy 
of the incident photon is close to the gap between the conduction and valence 
bands. We identify the diagram which gives the dominant contribution to Ra- 
man intensity in this regime and show that it can explain the unusual features 
in the two-magnon profile and in the two-magnon peak intensity dependence 
on the incoming photon frequency. 
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There is a widespread belief that strong electron-electron correlations in the high-T c com- 
pounds may hold a clue to the phenomenon of high-temperature superconductivity One 
of the manifestations of these correlations is the fact that the insulating parent compounds 
are antiferromagnets. An important probe of antiferromagnetism is magnetic Raman scat- 
tering [0,0. Its prominent signature in the underdoped high-T c materials is a strong peak 
observed at about 3000cm -1 . To first approximation, this peak can be attributed to inelastic 
scattering from the two-magnon excitations f||-|5|]. 

The traditional framework for understanding the two-magnon Raman scattering in anti- 
ferromagnets has been an effective Hamiltonian for the interaction of light with spin degrees 
of freedom known as the Loudon-Fleury Hamiltonian ||, H = a X)/y)(ei'Rij)(e/-Ry)Sj • Sj, 
where and e/ are the polarization vectors of the in- and outgoing photons, a is the cou- 
pling constant, and Rjj is the vector connecting two nearest neighbor sites % and j. Shastry 
and Shraiman [[?] have recently derived this Hamiltonian starting from the large-t/ Hubbard 
model. Working in a localized basis, they performed a hopping expansion controlled by 
tj (U — u), where t and U are the nearest-neighbor hopping and on-site Coulomb repulsion, 
and uj is of the order of the photon frequencies. The leading term in the expansion turned 
out to be the Loudon-Fleury Hamiltonian. 

This theory works well when the frequencies of the incoming and outgoing photons 
are considerably smaller than the gap between the conduction and valence bands, which 
is roughly 2eV. The experimental reality in high-T c materials is such, however, that the 
two-magnon scattering is measured mostly in the resonant regime, when the frequencies of 
the ingoing and/or outgoing photons are close to the gap value and the cross-sections vary 
strongly as the incident photon frequency is varied @,f|. Luckily, it is in this regime that 
the cross-sections sensitively depend not just on magnetic, but on the carrier properties as 
well, and this makes understanding the data particularly important. 

The profile of the Raman cross-section as a function of the transferred photon frequency 
and the behavior of the two-magnon peak height as a function of the incident photon 
frequency are shown in Fig. [j]. 
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The key experimental features that require explanation are: (a) the two-magnon peak is 
asymmetric, with the spectral weight shifted to higher frequencies; (b) the Loudon-Fleury 
Hamiltonian predicts no scattering in the A\ g configuration, whereas experimentally the 
resonant A\ g cross-section is about half of that in the B\ g geometry; (c) there is only one 
peak in Fig. |l]b — ordinarily one might expect two peaks, the so called ingoing and outgoing 
resonances HlQ|,[lT| ; (d) a comparison with the dielectric constant shows that the strength of 
the two-magnon Raman scattering in Fig. |l|b is at its maximum away from the band edge, 
in fact right at the upper end of the features in the optical data that can be interpreted as 
the particle-hole excitations between the lower and upper Hubbard bands. 

In this paper we develop a diagrammatic approach to Raman scattering valid in both 
nonresonant and resonant regimes. For small frequencies, our results are identical to those of 
Shastry and Shraiman. However, for \u — U\ ~ O(J), we find that the dominant contribution 
to Raman scattering comes from a diagram which is subleading in the nonresonant region. 
We will argue that this diagram accounts for most of the experimental features in Fig. [I]. 

We start with the one-band Hubbard model with H = — t Y^u,j)( c i,o- c j,<T + ^- c ) + 
U J2i n i-] n ii- I n the presence of the slowly varying vector potential A(x, t) the Hubbard 
Hamiltonian gets transformed to 

H = H A =°-^i q -A q + 0(A 2 ), (1) 

where j g = J2k §^ Lc l+ q /2a c fc-<?/2,o- is the current operator, and e^. = — 2t(cos k x + cosk y ) is 
the electron dispersion. The resonant part of the scattering matrix element Mr is obtained 
from the term linear in A in the 2nd order of perturbation theory [|IIJ]. In this process, 
a photon with energy Ui and momentum which can be safely set equal to zero, creates a 
virtual particle-hole state of the fermionic system which can emit or absorb two spin-waves 
with momenta k and — k before collapsing into an outgoing photon with the energy ujf (see 
Fig. |2|). Our primary goal is to calculate the dependence of this matrix element on the 
incident photon frequency. 



We use the spin density wave (SDW) formalism |12[] to describe the electronic state at 



half-filling and the excitations around it. In the SDW formalism one introduces a long- 
range order in S q = J2 k c \+ q ,a a a,p c k,p with q = Q = (71", 7r) and uses it to decouple the 
Hubbard interaction term. The diagonalization then yields two bands of electronic states 



(the conduction and valence bands) with E k = ±ye| + A 2 , where the 2A ~ U in the strong 
coupling limit that is assumed throughout this work. In terms of the conduction and valence 
band quasiparticle operators a ka and b ka , the current operator is interband to leading order 
in t/U, j q= o =>- J2k^t( a ka^ka + b\ a a ka ) . We also need the magnon-fermion interaction. 
Its derivation in the SDW formalism is straightforward, as the magnons are described as 



collective modes in the transverse spin channel ||12|| . The answer is, for S = 1/2 



#el-mag = X) ( a l,k a -a,k+ q el $aa 
k,q 

+ a^ k b„ a>k+q e\<& ab {k, q) + (a -> b) + H.c.) (2) 

where e q are the magnon operators, r] q ,rj q = ("S7 1 ) ^ ? 7g = ( cos Qx + cosq y )/2, and, to 
leading order in t/U, the vertex functions are given by 

®aa,bb{k, q) = ±(t k + tk+q)Vq + ( e fc ~ e k+q)rj q ] 

^abM( k ^) = 2A [VqTVq] ■ (3) 

In the situation when the photon frequencies are much smaller than A all the energy 
denominators are of order U, and the dominant diagrams for the Raman vertex are simply 
those with the largest numerators. It then follows from (|3|) that one has to consider processes 
in which the fermion changes bands while emitting a magnon. A representative diagram is 
shown in Fig. 0a. We collected all the leading order diagrams for the Raman matrix element 



and obtained 



M k = a H e i a e* fa [cosq a (A 2 , + fifj - 2A g /iJ , (4) 



a=x,y 



where v^/A? = fjq + Vqi V%, \ — Vq — Vqi an d a = 16t 2 A/(4A 2 — cu 2 ). This is exactly 
the expression which Shastry and Shraiman obtained in their derivation of the Loudon- 
Fleury vertex for the Hubbard model. Observe that within the Loudon-Fleury model, the 



scattering in the A ig geometry (e$ = = T^f-), vanishes, because in the A lg geometry 

the Loudon-Fleury and Heisenberg Hamiltonians commute with each other. On the contrary, 
in the B\ g scattering geometry (e* = e/ = the Raman vertex is finite, M^ 1 ~ 

+ \g)(cosq x — cosq y ) = cosg ^^ g -. The profile of the two-magnon Raman scattering in 
the Big geometry, including the final state interaction between magnons, has been studied 
several times in the literature HH. The two-magnon intensity has a narrow peak at the 
transferred frequency uj ~ 3J. 

A more careful treatment is, however, necessary in the resonant region, when the incom- 
ing photon frequency is close to the gap value, and one can no longer neglect the quasiparticle 
dispersion in the denominators. In fact, most of the experiments on Raman scattering have 
been performed in the frequency range where both Ui and Uf differ from U only to order J. 
In this situation, we found that the diagrams with the intraband fermion-magnon vertices 
become dominant, since they contain more resonant denominators. We analyzed these dia- 
grams and found that the most singular contribution to the Raman vertex comes from the 
one in Fig. |||b. The internal frequency integration in this diagram results in 

M (2) o-V-' (|k " (%* • g) [/Vfc-g ~ A g 6 fc ] 2 

M R -- 8 ^7^ otp , ak\(,, n w. iJ , ;x\t.,. n~W 
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(ui - 2E k + i8)(ui - Q q - E k - E k ^ q + i5)(u f - 1E k _ q + iS) ' 



where Vt q = 2 Jyl — 7^ is the semiclassical spin-wave frequency. A study of the integral 
shows that there is a region of and Uf where all three of the denominators in (^) vanish 
simultaneously, and the velocities v*. = dE k /dk and v k _ q are antiparallel to each other (oth- 
erwise, the integral over k vanishes). This phenomenon is known as a triple resonance |T0| , pT| . 
Via a combination of analytical and numerical techniques, we found that for relevant uji the 
triple resonance in the Raman vertex occurs only in a narrow range of the final photon 
energies Uf. The region of triple resonances is shaded in Fig. ^. 

It is important for our considerations that the triple resonance in the Raman vertex oc- 
curs only if both excited magnons are on the mass shell (only then is the second denominator 
in (^j) a half-sum of the other two). This is true only for the diagram which does not contain 
final-state magnon-magnon interactions. On the other hand, for S = 1/2, the dominant 



contribution to the conventional two-magnon peak at ~ 3J comes from the diagrams with 
magnon-magnon interactions ||. In this situation the Raman spectrum R(u) can be con- 
sidered as containing two independent peaks: one is due to the triple resonance in Mr in the 
shaded region in Fig. |3], which for most of the experimentally measured ajj is located close 
to 4J, and the other, at transferred frequency of about 3J, is due to the magnon-magnon 
scattering. Without considering in detail the effects of the fermionic damping, which smear 
the singularity in Mr, we cannot conclude which of the two peaks is stronger. The exper- 
iments indicate that the peak at 3J is stronger than that at 4J, and the enhancement of 
the Raman matrix element at larger transferred frequencies is responsible for the observed 
asymmetric shoulder-like behavior of the two-magnon profile. Suppose we now fix uj at the 
two-magnon peak frequency 3 J, as in Fig. p], and consider the variation of the peak ampli- 
tude as a function of the incident photon frequency uJi. Obviously, this amplitude will by 
itself have a maximum when the two peaks in R{uj) merge, i.e., when the u = const line 
intersects the region of triple resonances. From Fig. |^ we see that the intersection occurs in 
a very narrow region of uo,i close to u^ 1 ^ = 2A + 8 J, where the particle and hole are excited 
near the tops of their respective bands JTB|. We calculated Raman vertex in the vicinity of 



the intersection and found that it diverges (in the absence of damping) as 

Mr ~ 7 res i kfi ' ( 6 ) 

where uj\ es can be well approximated by u'- es = u] nax — (u — 2J) 2 /8J. The 3/2 power of the 
denominator in (§) is due to triple resonance, while the small factor in the numerator comes 
from the vanishing of the numerator in (|5]) right at the top of the band (i.e., at k = 0). In 
practice, the difference between uj\ es and iv™ 1 ^ can be neglected, and Eq. @ yields inverse 
square- root singularity in Mr, which implies a linear singularity in the Raman intensity, 
R~ \M R \ 2 ~ Iwi-wr"! -1 . 

Eq. (||) is a key result. In essence, we have found that the intensity of the two-magnon 
peak increases by an inverse linear law as one approaches the upper edge of the fermionic 
band. We emphasize that the singularity at the top of the band exists, due to the triple 
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resonance effects, despite the vanishing of the numerator in (|5|) at this point. 

We now discuss how these (and other) results are related to experiment. We listed the 
key experimental features in the beginning of the paper. Here we comment on each of them: 
(a) asymmetry of the two-magnon peak profile: our theory predicts that for u smaller than 
cj[ es the two-magnon peak profile should be asymmetric, with the "shoulder-like" behavior 
at frequencies close to u = 4J. This is consistent with experimental observations. In 
particular, the experimentally measured two-magnon profile in PrzCuO^ was analysed |14 



and found to contain two peaks, a two-magnon peak at 3000cm -1 , and a smaller one at 
4000cm -1 , which is precisely as expected from our calculations; (b) selection rules: the 
leading diagram in the resonance regime contributes to the scattering in both B\ g and A lg 
geometries. The signals in both geometries have been observed in the experiments. Recall 
that the Loudon-Fleury theory predicts scattering only in the B\ g geometry; (c) a single 
peak: our theory predicts a single maximum in the two-magnon peak intensity measured as 
a function of the incident photon frequency, while from the "Loudon-Fleury" diagrams we 
might have expected two peaks, one at Ui = 2 A, and the other at ujf = 2 A (the incoming and 
outgoing resonances) |TIJ; (d) peak location and shape: our theory predicts the maximum 
of the two-magnon peak intensity measured as a function of u>i right near the upper edge of 
the quasiparticle fermionic band. This is consistent with the measurements of the dielectric 
constant, which show that the Raman scattering is strongest right at the upper edge of those 
features in the optical data that can be interpreted as particle-hole excitations between the 
lower and upper fermionic bands. We fitted the data on the peak intensity from Fig. [I] by 
our Eq. (^) and found a satisfactory agreement with the predicted inverse linear dependence, 
except in the immediate vicinity of the resonance, where the effects of fermionic damping 
become relevant. 

To summarize, we developed a diagrammatic approach to Raman scattering in antifer- 
romagnetic insulators which can be used in both the resonant and nonresonant regimes. 
We described for the first time the two-magnon Raman scattering in the resonant regime, 
when the incident and final photon frequencies are only O(J), apart from the gap between 



conduction and valence bands. This frequency range is relevant to recent experiments on 
undoped high-T c compounds. We identified the diagram which gives a dominant contribu- 
tion to the Raman vertex in this regime, and found the region in the (ui,Uf) plane where 
the Raman vertex is strongly enhanced due to triple resonance. We demonstrated that the 
triple resonance, combined with the SDW dispersion relation for the carriers, explains the 
unusual experimental features in the two-magnon profile and in the two-magnon peak in- 
tensity dependence on the incoming photon frequency. In particular, our theory predicts 
the maximum of the two-magnon peak intensity right at the upper edge of the features in 
the optical data, as observed in several materials 0. This serves as a partial verification 
of the SDW picture for the carriers, which, despite much theoretical work, has not been 
well-established experimentally in these materials. 

Beyond the scope of the present theory are the unexpectedly large width of the symmetric 
part of the two-magnon peak, which is probably related to the magnon damping due to the 
interaction with phonons ||15| , and the existence of a considerable Raman signal R(u) well 
above the maximum possible two-magnon energy (i.e., 4 J) |16|| , which may be related to 
chiral spin fluctuations |T7| . 
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FIGURES 

FIG. 1. (a) A typical Raman cross-section in Y B<iCu^O% as a function of transferred photon 
frequency. A two-magnon peak is clearly seen, (b) The strength of the two-magnon peak as a 
function of incoming photon frequency. Also shown is the imaginary part of the dielectric constant. 
Data courtesy of the authors of Ref . || . 

FIG. 2. (a) A representative diagram which contributes to the Loudon-Fleury Hamiltonian 
at small incident frequencies. Each fermion can belong to either the valence (dashed line) or 
conduction (solid line) band. The emitted magnons are denoted by the solid wavy lines, (b) The 
most singular diagram at resonance. 

FIG. 3. The triple resonance region (shaded) in the (uJi,to) plane where oj = Ui — ut. The 
horizontal line corresponds to the position of the two-magnon peak which for definiteness we chose 
to be at u) = 2.8 J which is the value one obtains in the 1/S expansion neglecting the renormalization 
of J |J]. 
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